AD-AO03  173 
UNCLASSIFIED 


STANFORD  UNIV  CA  DIV  OF  APPLIED  MECHANICS  F/G  11/4 

ELASTICITY  THEORY  OF  COMPOSITES. (U) 

MAR  80  J  R  WILLIS  N000U-76-C-0054 

SUDAM-80-2  NL 


Office  of  Naval  Research 
Department  of  the  Navy 
Scientific  Report 
Contract  N00014-76-C-0054 


72/ 


DIVISION 


APPLIED 

MECHANICS 


DEPARTMENT 


ELASTICITY  THEORY  OF  COMPOSITES 


MECHANICAL 

ENGINEERING 


J.  R.  Willis 


I  i 


it-  .V 

i 


SUDAM  Report  No.  80-2 


March  1980 


This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


STANFORD 

UNIVERSITY 


STANFORD. 

CALIFORNIA 

94305 


DTIC 

k ELECT  I 


APR  1  8  1980 


80  4  17  0  16 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (TWin  D«l«  gnfuQ 


DO  1  jan*7S  1473  EDITION  OF  I  NOV  88  IS  OBSOLETE 

/  /s  y  .  9 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  page  (■IWifu  D  «•  Inland) 


ELASTICITY  THEORY  OF  COMPOSITES 


by 

J.  R.  Willis 


SUBJECT  INDEX 


Correlation  functions  33 
Cracked  solid  61-62 
Fibre-reinforced  media 

long  fibres  55-58 
short  fibres  58-61 
Finite  deformations  9 
Green's  function 

dynamic  65,71 
static  10,20 
Hierarchy  equations  68 
Inclusions 

general  equation  21 
solution  for  ellipsoid  21-22 
solution  for  elliptic  cylinder  24-25 

Momentum  polarization 
definition  65 
operator  equation  for  66 
Overall  compliances  42-43 
Overall  moduli 

basic  definitions  4-6 
bounds  of  third  order  43-46 
estimate  to  order  C*  46 
Hashin-Shtrikman  bounds  41 
of  dilute  suspension  25 
Reuss  estimate  8 

self-consistent  estimates  27-31,  47-48 
symmetry  of  7 
Voigt  estimate  7 
Plane  waves  70 
Polycrystal  52-55 
Porous  medium  50 

Quasicrystalline  approximation  38,69 
Radon  transform  21 
Random  media  32-36 
Renormalization  63 
Smoothing,  method  of  63 
Strain  polarization  16 
Stress  polarization 
definition  12 

operator  equation,  displacement  boundary  condition 
operator  equation,  dynamic  problems  66 
Variational  principle 

complementary  energy  7 

for  dynamic  problems  67 

Hashin-Shtrikman  13-14 
minimum  energy  8 


ELASTICITY  THEORY  OF  COMPOSITES 


by 

J.  R.  Willis 
Bath  University 


SUMMARY 

The  determination  of  the  overall  elastic  moduli  of  random  composites 
is  approached  using  variational  principles.  Information  on  the  configu¬ 
ration  of  a  composite  is  never  complete  and  strict  estimates  can  do  no 
more  than  take  the  form  of  bounds  that  use  whatever  statistical  inform¬ 
ation  is  available.  The  classical  Hashin-Shtrikman  bounds  are  deduced 
as  special  cases  of  bounds  that  allow  explicitly  for  any  two-point  cor¬ 
relations  and  it  is  shown  that  these  general  bounds  take  a  simple  form 
not  given  before  for  any  two-phase  composite.  Self-consistent  estimates 
are  related  to  the  variational  formulation.  Bounds  of  higher  order  are 
discussed  and  equations  that  define  a  new  optimal  bound  of  third  order 
are  given  explicitly.  The  concluding  sections  treat  wave  propagation 
problems  from  the  standpoint  of  a  new  variational  principle  which  reduces 
to  the  Hashin-Shtrikman  principle  in  the  static  limit. 
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1.  INTRODUCTION 

The  subject  of  this  article  is  the  theoretical  determination  of  the 
overall  properties  of  inhomogeneous  materials  such  as  fibre-reinforced 
plastics.  The  distinguishing  feature  of  these  materials,  to  which  the 
generic  term  'composite'  is  applied,  is  that  they  are  strongly  inhomo¬ 
geneous  relative  to  a  small  length  scale  (the  microscale)  and  yet,  rela¬ 
tive  to  a  larger  scale  (the  macroscale)  over  which  variations  in  applied 
loads  are  significant,  they  appear  either  homogeneous  or  else  their 
'average'  properties  appear  to  vary  smoothly.  Generally,  a  composite 
may  comprise  n  different  phases,  distributed  either  deterministically  or 
at  random,  and  it  will  be  assumed  that  the  microscale  defined  by  the 
local  phase  geometry  is  sufficiently  large  for  each  phase  to  be  treated 
as  a  continuum:  this  is  certainly  true  of  fibre-reinforced  plastics  and 
it  is  true  also  of  polycrystalline  aggregates.  The  latter  represent  a 
limiting  case  in  which  the  number  of  phases  is  indefinitely  large,  each 
crystal  orientation  being  regarded  as  defining  a  different  phase. 

The  article  will  be  concerned  with  the  determination  of  the  overall 
elastic  response  of  a  composite,  since  this  particular  aspect  of  material 
behaviour  can  be  treated  with  reasonable  precision.  Some  of  the  methods 
to  be  described  have  wider  application  (when  suitably  generalized) ,  for 
example  to  the  elastoplastic  response  of  a  polycrystal,  but  study  of  the 
simpler  elastic  problem  permits  some  assessment  of  limitations  inherent 
in  these  methods,  by  comparison  with  results  obtained  from  methods  that 
allow  greater  precision  but  are  less  easily  generalized. 

The  elastic  properties  of  composites  have  been  studied  quite  inten¬ 
sively  during  the  past  twenty  years  and  many  of  the  simpler  problems  are 


at  a  definitive  state  of  development.  A  reasonably  comprehensive  review 
is  contained  in  the  recent  book  by  Christensen  (1979).  The  major  gaps 
in  Christensen's  presentation  concern  allowance  for  random  phase  geo¬ 
metry  when  this  is  not  isotropic,  and  wave  propagation  in  randomly  inhomo¬ 
geneous  composites.  Progress  has,  in  fact,  been  made  in  both  of  these 
areas  which  happen,  too,  to  be  of  particular  interest  to  the  writer. 
Correspondingly,  attention  is  focussed  in  this  article  upon  problems  for 
random  media,  which  need  not  be  isotropic.  The  presentation  is  self- 
contained;  the  longer-established  results  are  deducible  as  special  cases 
and  some  of  the  more  important  of  these  are  discussed  explicitly,  though 
not  in  the  detail  given  by  Christensen.  Apart  from  introductory  material 
the  present  review  is  largely  complementary  to  one  recently  prepared 
by  the  writer  (Willis,  1980d) .  This  covered  a  somewhat  broader  range 
of  topics  and  emphasised  more  strongly  composites  consisting  of  a  matrix 
reinforced  by  inclusions.  Distinctive  features  of  the  present  work  are 
the  discussion  of  bounds,  including  a  new  simple  formula  which  allows 
for  arbitrary  two-point  correlations  in  a  two-phase  composite,  new 
optimal  bounds  of  third  order  and  a  variational  treatment  of  wave  prop¬ 


agation  problems. 
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2.  OVERALL  PROPERTIES 

Underlying  the  'overall  modulus'  concept  is  the  belief  that,  in  a 
specimen  subjected  to  some  system  of  loading,  the  stress  and  strain  fields 
vary  in  a  complicated  fashion  from  point  to  point  but,  if  some  form  of 
'running  average'  is  taken  which,  in  effect,  replaces  the  fields  at  X 
by  weighted  means  evaluated  over  some  neighbourhood  centred  on  x  , 

•AT 

whose  dimensions  are  small  relative  to  overall  specimen  dimensions  and 
yet  large  relative  to  the  microscale,  then  the  'averaged'  stress  and  strain 
fields  will  vary  smoothly  relative  to  the  macroscale.  It  is  difficult 
to  define  the  averaging  procedure  precisely  for  a  general  composite  but 
it  corresponds,  hopefully,  to  what  might  be  measured  by  a  strain  guage 
or  a  transducer.  The  details  of  the  fields  in  the  neighbourhood  of  x 
are  then  expected  to  be  determined  just  by  their  mean  values  at  x*  and 
the  local  configuration  of  the  material.  These  expectations  are  indeed 
borne  out  for  composites  with  a  periodic  structure,  for  which  the  stress 
and  strain  fields  may  be  expressed,  asymptotically  at  least,  as  products 
of  periodic  functions  with  the  period  of  the  structure,  and  amplitudes 
that  vary  smoothly  on  the  macroscale.  For  static  problems,  these  ideas 
were  enunciated  by  Sanchez  Palencia  (1974).  The  Floquet  theory  approach 
of  Kohn,  Krumhansl  and  Lee  (1972)  for  dispersion  relations  for  periodic 
composites  is  essentially  similar,  except  that  it  does  not  even  require 
the  wavelength  associated  with  the  'amplitude'  function  to  be  large. 

Overall  moduli  are  defined  so  that  they  relate  'average'  stress  to 
'average'  strain.  Symbolically,  if  stress  and  strain  tensors  are  denoted 
•  t  and  their  local  averages  are  correspondingly  represented  as  •*  ,  8  , 
then  the  tensor  of  overall  moduli  L  is  defined  so  that 
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=  L  e  . 


(2.1) 


Strictly,  if  the  relation  (2.1)  is  to  apply  when  ,  ft  are  generated 
from  arbitrary  boundary  conditions,  then  L  should  be  a  non-local  oper- 
ator  (Kroner,  1977)  but,  asymptotically,  when  boundary  conditions  are 
such  that  ft*  ,  ft  vary  smoothly  relative  to  the  microscale  (and  so  hardly 
at  all  in  a  neighbourhood  of  *  ) ,  L  reduces  to  a  tensor  of  moduli  ex- 
cept,  of  course,  in  a  thin  layer  close  to  the  boundary  of  the  specimen 
which  is  excluded  from  detailed  consideration. 

If  the  validity  of  (2.1)  is  accepted,  the  overall  moduli  L  may  be 
determined  experimentally  by  studying  the  response  of  a  macroscopically 
uniform  specimen  when  it  is  subjected  to  loads  that  generate  uniform 
average  stress  and  strain  fields.  In  fact,  the  tensor  of  moduli  L.  may 
be  defined  in  an  unambiguous  manner  by  subjecting  a  specimen  that  occupies 
a  region  \/  with  boundary  to  surface  displacements  that  have  the 

form 


U'  s  JL, 

1  j)  J 


(2.2) 


relative  to  a  Cartesian  basis,  the  summation  convention  being  implied  for 
repeated  suffixes.  Application  of  Gauss'  theorem  then  shows  that  the 
volume  average  of  the  strain  actually  produced  throughout  ^  is  C 

precisely  (ft.  having  components  ft..  ).  If  the  average  over  V  of  the 

~  \ 

stress  C*  is  denoted  by  ft",  then  depends  linearly  upon  ft  and  de- 

*•»  __«*<•  ■*» 

fines  the  tensor  L.  in  (2.1).  The  above  discussion  assumes,  of  course, 


that  the  phases  of  the  composite  are  linearly  elastic  and  that  there  is 
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perfect  adhesion  across  phase  boundaries.  It  may  be  noted  that  the  com¬ 
ponents  •TJj  of  the  mean  stress  may  be  determined  from  measured  values 
of  surface  traction,  through  the  relation 


(2.3) 


where  \/  has  also  been  used  to  denote  the  volume  of  the  region  ^  to 
avoid  further  symbols,  and  /?.  ,  the  normal  to  having  compo- 

J  J*  * 


nents  rt. 


Equation  (2.3)  follows  from  Gauss'  theorem,  together  with 


the  requirement  that  the  stress 


should  be  self-equilibrated;  that  is. 


=  O  (2.4) 

except  at  interphase  boundaries,  across  which  ft.  is  continuous. 

J 

The  latter  proviso  is  included  in  (2.4)  if  the  derivatives  are  inter¬ 
preted  in  the  sense  of  generalized  functions  (Gel'fand  and  Shilov,  1964): 
this  interpretation  will  be  used  throughout  the  sequel  to  avoid  the  need 
for  explicit  discussion  of  jump  conditions  and  it  should  not  be  assumed 
that  sharp  changes  in  properties  are  thereby  excluded. 

An  alternative  definition  for  L  is  obtained  by  applying  to  the 

•»>» 

specimen  the  traction  boundary  condition 


(2.5) 


Equation  (2.3)  then  ensures  that  the  average  of  the  actual  stress  over 


y  is  and  L  follows  from  (2.1),  this  time  with  &  calculated  as 
the  average  of  the  strain  over  V  .  Again,  can  be  determined 
from  surface  values  of  displacement  alone,  by  use  of  the  identity 


=  _U  [  (w- w,  ♦  u,H:)  US, 


(2.6) 


Thus,  either  experiment  determines  from  values  that  are  accessible 
to  measurement,  in  principle  at  least.  They  were  both  proposed,  as  pre¬ 
cise  mathematical  definitions  for  L  ,  by  Hill  (1963a).  The  two  defini- 
tions  are  not  equivalent  to  one  another  for  arbitrarily  inhomogeneous 
bodies  and  the  extent  to  which  they  agree  in  their  estimates  for  (_  pro- 
vides  a  partial  check  on  the  validity  of  the  ’overall  modulus'  concept. 

The  definitions  are  interesting  for  another  reason,  also  developed 
by  Hill  (1963a),  that,  when  either  of  the  boundary  conditions  (2.2)  or 
(2.5)  is  applied,  then  the  average  energy  density  in  the  specimen  is 
given  by 


2  U  s  <s-  e  =  e  L 


(2.7) 


precisely,  again  by  application  of  Gauss'  theorem  combined  with  one  of 
(2.3),  (2.6).  The  shorthand  notation  implicit  in  (2.1)  has  been  extended 
in  (2.7)  which,  in  suffixes,  would  read 


V 


(2.8) 


In  the  sequel,  suffixes  will  be  avoided  wherever  possible. 
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The  derivation  of  (2.7)  has  assumed  the  usual  symmetry  of  the  tensor 

of  moduli  L  which  applies  to  the  actual  inhomogeneous  specimen.  Cor- 

responding  symmetry  of  1.  may  be  demonstrated  by  calculating  the  incre- 
- 

ment  su  in  mean  energy  corresponding  to  an  increment  SI  in  mean 
strain:  the  argument  that  led  to  (2.7)  gives 


$U  =•  k  & 


(2-9) 


which,  when  compared  with 


the  desired  result. 


calculated  by  differentiating  (2.7),  gives 


Bounds  upon  components  of  the  tensor  of  overall  moduli 


may  now 


be  developed  from  (2.7)  and  the  classical  variational  principles  for 
elasticity.  For  the  boundary  condition  (2.2),  the  minimum  energy  prin¬ 


ciple  asserts  that 


ZU  < 


if  t*  L 

vl. 


e*  dx 


(2.10) 


for  a  strain  field  er  derived  from  any  displacement  field  14  that 

conforms  with  (2.2)  on  w .  The  simplest  field  that  can  be  substituted 

jf  „ 

into  (2.10)  is  "C-  s  'Gs  ;  this  gives 


k  -  Lv ) 


(2.11) 


where  Ly  denotes  the  average  of  L.  over  }/  ;  it  was  introduced  by 
Voigt  (1889).  The  symbol  *  in  (2.11)  defines  an  ordering  in  the 

sense  that  the  quadratic  form  £.  I  L  L  is  positive  semidefinite 

*  *  «»r  «w  /•* 


is  positive  semidefinite. 


for  symmetric  second-order  tensors  .  Dually,  for  the  boundary  condi- 

•v* 

tlon  (2.2),  the  complementary  energy  principle  asserts  that 


2-f  b*  U.  dS  -  L  [ cr* n  <r*  Mx  <  Z  U 

V  L,  ~  ~  vJ, - - 


(2.12) 


'W  V 


,  s  H;  ,  tf"  is  any  self-equilibrated  stress  field, 

ji  *J  J  ~ 

U  =  —  'C. ..  Jt-  and  M  is  the  tensor  of  compliances,  inverse  to  L.  . 

*  J  —  ~ 

The  simplest  stress  field  C5"  to  substitute  into  (2.12)  is  a  constant; 
this  given,  upon  use  of  Gauss'  theorem. 


<■  2U, 


(2.13) 


where  is  the  average  of  M  over  \/  ,  introduced  by  Reuss  (1929). 

««%» F\  /Vv 

The  left  side  of  (2.13)  is  maximized  by  taking 


= 


(2.14) 


where  L  is  inverse  to  M M  •  It  follows  then  that 


<  L  <  Ly . 


(2.15) 


The  result  (2.15)  is  due  to  Hill  (1952).  It  can  also  be  deduced  in  a 
corresponding  fashion  if  L  is  defined  by  the  boundary  condition  (2.5). 

«w 

The  bounds  L  ,  Lw  depend  only  upon  the  moduli  and  concentrations  of 
the  phases  and  take  no  account  of  how  the  phases  are  distributed.  The 
construction  of  better  bounds  which  make  more  allowance  for  the  structure 
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of  the  composite  is  considered  in  later  sections.  Before  proceeding, 
however,  it  may  be  remarked  that  corresponding  mean  value  theorems  can 
be  developed  for  a  body  subjected  to  finite  deformation.  The  relation 
(2.3)  is  valid  in  any  continuum  but  it  is  useful  to  have  a  description 
that  employs  Lagrangian  variables,  based  upon  a  fixed  reference  state. 
Relative  to  such  a  state,  Hill  (1972)  has  shown  that  an  equation  of  the 
form  (2.3)  applies  if  the  Cauchy  stress  ©”  is  replaced  by  nominal  stress 

and  X.  is  interpreted  as  a  Lagrangian  coordinate.  Alternatively, 

•1 

may  be  replaced  by  Kirchhoff  stress  and  X L  is  interpreted  as  the  cur¬ 
rent  (Eulerian)  coordinate,  expressed  as  a  function  of  the  Lagrangian 
coordinates.  An  equation  like  (2.6)  (not  symmetrized)  applies  to  the 
deformation  gradient,  and  mean  Kirchhoff  stress,  nominal  stress  and  de¬ 
formation  gradient  are  related  exactly  as  though  the  fields  were  uniform. 
Similar  relations  apply  to  rates  which  allow,  for  elastic  materials,  gen¬ 
eralizations  of  (2.7)  and  (2.9);  some  implications  for  nonlinearly  elas¬ 
tic  composites  have  been  explored  by  Ogden  (1978) . 
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3.  STRESS  POLARIZATION  AND  VARIATIONAL  PRINCIPLES 

Bounds  better  than  (2.15)  require  the  construction  of  stress  and 
strain  fields  that  make  some  explicit  allowance  for  the  microscopic 
arrangement  of  the  composite.  For  this  purpose,  it  proves  useful  to 
introduce  a  homogeneous  'comparison'  material  with  moduli  L.  and  to  set 


L  =  L,  *fL . 


(3.1) 


Substitution  of  (3.1)  into  the  constitutive  relation 


=  L 


(3.2) 


and  thence  into  the  equilibrium  equation  (2.4)  gives 


aiv(L.c)  =  o  .  xeV, 


(3.3) 


which  is  to  be  solved  together  with  whatever  boundary  condition  is  given. 
The  discussion  to  follow  will  be  phrased  in  terms  of  the  displacement 
condition 


14.  S  u  X  £  ~)V 

Z.  £©  >  ±  i 


(3.4) 


with  a  mind  to  subsequent  specialization  to  (2.2).  If  the  Green's  func¬ 
tion  Q  for  the  homogeneous  comparison  material  is  known,  (3.3)  and 
(3.4)  may  be  expressed  together  in  the  form 


..  . _ ... 
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u 


(SL  ft  )(x')  Mx't 


(3.5) 


where  s  the  field  that  the  boundary  condition  (3.4)  would  gen¬ 

erate  by  itself  in  the  homogeneous  comparison  body.  An  integral  equation 
for  the  strain  t(z)  is  now  deduced  by  integration  by  parts,  followed 
by  differentiating  with  respect  to  X  .  This  gives 


V 

where  C>  „  is  the  strain  associated  with  14  and  I  (x,t/) 

f7  (xy)  =  ~>l (* , *')/)*< I 

tJK1-1-7  /*  ’  tlf.. 


(3.6) 


has  components 


(3.7) 


It  has  a  singularity  of  order  which  is  to  be  interpreted 
in  the  sense  of  generalized  functions  (Gel ' fand  and  Shilov,  1964).  In 
the  sequel,  (3.6)  will  be  written  even  more  briefly  as 


■e 


(3.8) 


If  the  components  of  0  L  are  small,  equation  (3.8)  can  be  solved  by 

M 

iteration  and  the  resulting  series,  suitably  truncated,  can  be  substituted 
into  the  energy  and  complementary  energy  principles.  This  has  been  de¬ 
veloped  in  detail  by  Dederichs  and  Zeller  (1973);  the  validity  of  the 
bounds  is,  of  course,  independent  of  convergence  of  the  perturbation 
series. 

It  is  possible,  however,  to  make  rather  better  use  of  the  work 
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implicit  in  the  construction  of  the  bounds  of  Dederichs  and  Zeller  (1973), 
by  considering  a  form  for  that  contains  some  parameters  which  may  be 
chosen  optimally.  It  is  helpful  for  this  purpose  (as  well  as  for  other, 
subsequent,  developments)  to  introduce  the  stress  polarization  "V  so 
that 


X  ■  -  (k-!=o)t 


(3.9) 


and 


<5~  =  L  e  +  r  . 


(3.10) 


Equation  (3.8)  gives 


e  * 


Ft 


(3.11) 


and 


,  for  any  choice  of  TT  ,  equations  (3.11)  and  (3.10)  generate 


strain  and  stress  fields  6-  and  that  are,  respectively,  compatible 

with  boundary  conditions  and  self-equilibrated.  If  the  mean  energy  den¬ 
sity  in  the  actual  body  is  called  W  ,  the  minimum  energy  principle  gives 


2U  *&*, L.s)  +(£,*€), 


(3.12) 


where  the  'inner  product'  is  defined  so  that 


(3.13) 
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The  first  term  on  the  right  side  of  (3.12)  is  simplified  by  employing 
(3.10)  and  (3.11)  for  •“*,  t?  and  IT*  and  using  the  'virtual  work' 
equality 
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This  principle  was  deduced  from  the  field  equations  by  Hashin  and  Shtrikman 
(1962a)  and  directly  from  (3.16)  by  Willis  (1977).  Here,  however,  it 
seems  particularly  apt  to  give  a  derivation  which  corresponds  to  that 
of  Hill  (1963b).  This  shows  how  the  Hashin-Shtr ikman  principle  is  re¬ 
lated  to  the  classical  energy  principles.  Reverting  to  (3.15),  there¬ 
fore,  the  offending  last  term  is  rewritten 

(rvsk  Cl")  -■  ('[csij,*fJr-  (GH)"+rJr'-2*)) 


since  the  exact  T  satisfies  (3.16),  the  first  term  on  the  right  of 
(3.18)  is  zero  when  31  anc*  *s  smaH  when  ^  ^  is  a  good  approxi¬ 

mation  to  £  .  Hence,  (3.15)  gives 


19) 


where 


£*  has 


been  written  for  the  "error’ 


RsLyVrjr-  *. 


If,  now,  i-s  chosen  in  such  a  way  that  SL  is  negative  definit 


re 


e  at 


each  point  of  V  ,  discarding  the  last  term  in  (3.19)  gives 


a(U-U.)  ^  -  Rl*). 


(3.20) 


Hi")  w  Wo~  U  )  ,  whatever  L  ^  is  chosen. 


The  extreme  value  of 
The  Inequality  (3.19)  shows  that  the  Hashin-Shtr ikman  principle  embodied 
in  (3.20)  is  weaker  titan  the  corresponding  classical  principle,  through 
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neglect  of  the  terra  quadratic  in  C *  .  Hill  (1963b)  made  no  use  of  the 
operator  equation  (3.16)  but  instead  gave  in  the  form 


e*.  (SLy'te*-  Ltf). 


The  complementary  energy  principle 

z(z*  £.)-  (T*  Mf")  £  2W 


may  be  expanded  similarly:  with  the  notation 


M  *  M„  ♦  SM 


so  that  is  inverse  to  L  . ,  and  some  use  of  the  identity 

a.  <*  *• 

L#*Q  -  fM,*  S-M)  - 

it  follows  that 


(3.23) 


(3.24) 


(3.25) 


If  (_  is  chosen  so  that  SM  is  negative  definite  (or,  equivalently, 
so  that  5L  is  positive  definite),  (3.25)  implies 


-Jit*)  <Z  z(U-U.). 


The  traction  boundary  value  problem  can  be  treated  in  a  similar  way. 


Its  duality  with  the  displacement  problem  is  best  displayed  by  intro¬ 
ducing  a  strain  polarization  ,  so  that 


%  =  g.i  t  =  a,*:  -  ? 


(3.27) 


If  tractions  are  specified  that  would  generate  the  stress  field  jT’g  in 

the  comparison  material,  an  operator  A.  may  be  defined  so  that 

I 


6“  =  <r  +  A  7 


(3.28) 


which  is  analogous  to  (3.11).  The  corresponding  analogue  of  (3.16)  is 

[4,  +(X!1)',J?  "-fo,  (3-29> 

with  which  the  functional 

fM'Hl'UsWh')  -  2Cf>  r.) 

may  be  associated.  The  minimum  energy  principle  and  the  complementary 
energy  principle  now  in  effect  change  roles.  The  analogue  of  (3.12), 
from  the  complementary  principle,  is 

U 

which,  expanded  in  terms  of  *1  ,  is 
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(3.32) 


-  U.)  £-(  f,  d,  f)  *  {<£„  S  M  r.)  ♦  a  fe.,  ) 

This,  finally,  can  be  subjected  to  the  reasoning  that  produced  (3.19), 
to  give 


(3.33) 


so  that 


*  'ft!*) 


(3.34) 


whenever  M  is  chosen  so  that  SM  is  negative  definite  (or,  equiva- 
lently,  so  that  ^L.  is  positive  definite).  The  "error"  £*  in  (3.33) 
still  conforms  to  (3.21);  in  terms  of  y  ^  , 


t.l*  *-fVWb'-r 


o  • 


(3.35) 


It  is,  of  course,  possible  to  express  the  functional  in  terms 

of  and  in  terms  of  ,  with  ^  and  related  by  (3.27)^. 

For  example,  equation  (3.16)  can  be  written 


[UQ"L  -  t.CLo]? 


(3.36) 


from  which  follows,  using  a  rearrangement  of  (3.24),  an  equation  of  the 
form  of  (3.29),  with  the  operator  ^  replaced  by  ^  ,  where 
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a  «  l  -  l  r  L  . 

Mo  M,® 


(3.37) 


by 


Correspondingly ,  fa*) -pi*)  ,  where  m  is  defined 

but  with  replaced  by  A  •  Dually,  where 

-*  p  »  *  n 


(3.30) 


is  defined  by  (3.17),  with  f  replaced  by  ,  where 


r  =  H  -  M  AM 

I  0  '*-1 


(3.38) 


Finally,  there  are  duals  of  (3.25)  and  (3.26),  obtained  from  the  minimum 
energy  principle  for  the  traction  problem: 


(3.39) 


and 


TL 


whenever  $  L  is  negative  definite  everywhere  in 


in  V 


(3.40) 
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4.  INCLUSION  PROBLEMS 

For  any  composite  that  comprises  a  matrix  reinforced  (or  weakened) 
by  inclusions,  a  fundamental  problem  is  to  determine  the  way  in  which 
any  particular  inclusion  perturbs  the  field  that  the  matrix  would  support 
in  its  absence.  Interaction  effects  between  inclusions  are  significant 
but  the  problem  of  finding  the  effect  of  a  single  inclusion  in  an  infi¬ 
nite  homogeneous  matrix  provides  an  essential  building  block  for  consider¬ 
ing  a  general  composite.  This  basic  problem  is  discussed  now.  The  inte¬ 
gral  equations  (3.8)  or  (3.16)  are  relevant  and  it  is  natural  to  identify 


the  comparison  material  with  the  matrix.  With  a  slight  distortion  of 
earlier  notation,  therefore,  the  matrix  is  taken  to  have  moduli  L & 
while  the  inclusion,  which  occupies  a  region  J)  ,  is  taken  to  be  uniform, 
with  moduli  L  .  The  field  ft  -  is  now  the  strain  field  that  would  be 

a. 

present  if  the  matrix  contained  no  inclusion  and  the  term  involving  f* 

Mm 

represents  the  perturbation  of  this  field  produced  by  the  inclusion.  SL 
is  zero  except  over  T)  and  (3.8)  is  now  an  integral  equation  only  for 

XCUi  elsewhere,  it  provides  a  representation  for  the  strain  field  . 

•%» 

Equivalently,  is  non-zero  only  over  2  and  equation  (3.16)  applies 

only  for  X  £  P  .  Either  of  these  equations  could  be  studied  equally 


well  but  there  is  some  later  advantage  in  concentrating  upon  (3.16). 

The  region  "P  occupied  by  the  inclusion  will  be  assumed  to  be 
bounded.  The  three-dimensional  problem  will  be  considered  first  but 
afterwards  a  short  discussion  will  be  given  of  the  corresponding  two- 
dimensional  problem  for  which  the  inclusion  may  be  an  infinite  cylinder. 
The  matrix  is  infinite  in  extent  and  so  the  operator  f  that  appears  in 
(3.16)  is  obtained,  via  (3.7),  from  the  infinite-body  Green's  function 
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C |  .  This  satisfies,  in  matrix  notation,  the  differential  equation 


lJv)^*  I =  £  , 

MV  ^ 


(4.1) 


where 


LW 


is  a  differential  operator  whose  symbol 


L  W 


components 


(4.2) 


Equation  (4.1)  can  be  solved  by  noting  the  plane-wave  decomposition 


w.j.-j-J  «'Ui) 


(4.3) 


given  by  Gel'fand  and  Shilov  (1964).  Solving  the  differential  equation 


(4.1)  with 


U*) 


replaced  by 


S  "(}.*) 


and  then  performing  the  super¬ 


position  implied  by  (4.3)  gives 


-S  [LJDl'Vd.  x) 


(4.4) 


l)|*l 

The  components  of  (so  called  to  correspond  with  (7  )  now  follow 

using  (3.7): 

w>J 45  UU)£i,cn  j«» 

1  ”  |)  |=,  S)(M) 


and  substitution  of  (4.5)  into  (3.16)  gives 
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n  JM|*i 


51  £<,(*),  £  *  <4'6) 


where 


i V i)  *  *  i)  - 1  \tj*') 


(4.7) 


^  ca-^ec*  t^le  transform  of  T  In  (4.6),  the  term 

represents  the  part  of  the  integrand  of  (4.5)  that  depends  upc 
J  only;  it  is  actually  the  Fourier  transform  of  /""* 


It  is  not  clear  whether  (4.6)  is  superior  to  (3.16)  in  general  but 


it  is  particularly  useful  when  J)  is  an  ellipsoid.  To  illustrate  this, 
consider  first  the  case  for  which  J)  represents  a  sphere  of  radius  ft  , 
centred  at  the  origin.  If  T  is  taken  constant  over  ^  ,  its  Radon 
transform  is  obtained  as  times  the  area  of  intersection  of  the  plane 
!,  Xs  j>  with  J)  ;  that  is, 


(4.8) 


The  second  derivative  that  appears  in  (4.6)  is  then  constant  when  |5  |s  I 
and  1C  £  J)  >  so  that  the  left  side  of  (4.6)  is  constant  over  T)  .  If 

is  taken  constant,  therefore,  t  may  be  taken  constant  and  must 


satisfy  thu  algebraic  equation 


[(HT'*£]t  =  i.. 


(4.9) 


where  the  constant  tensor  I  is  given  by 


r 

I 


l: 
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P  =  —  f  *s  p^j ) 

«r  1 

This  result  was  first  given  by  Kneer  (1965), 
the  ellipsoid 

P  *  J  x  :  xT  A^Ax  <  J  J 


(4.10) 


The  case  for  which  Is  is 


(4.11) 


can  be  reduced  to  the  one  just  considered  by  the  dual  transformations 


1  ■-  A  *  ,  5  -  AT  S  ; 


(4.12) 


this  reproduces  equation  (4.9)  except  that 


is  now  given  as 


P  =  ± 


4n\h\ 


f  <sf 7£)[ 

m=i 


(4.13) 


This  result  was  given  first  by  Khatchaturyan  (1967)  for  a  generally 
anisotropic  matrix.  For  an  isotropic  matrix,  the  result  was  derived 
using  methods  of  potential  theory  by  Eshelby  (1957),  who  expressed  P 
in  terms  of  elliptic  integrals.  Eshelby  (1961)  showed  further  that,  for 
an  inclusion  in  an  isotropic  matrix,  a  polynomial  field  6-  ( X  )  would 

generate  a  polynomial  field  within  the  inclusion:  the  corresponding  re¬ 
sult  for  a  generally  anisotropic  matrix  was  given  by  Asaro  and  Barnett 
(1975)  and  Willis  (1975),  essentially  by  demonstrating  that,  if  T  j 
is  a  polynomial  of  degree  O  in  X  ,  then  T  }  j  is  a  polynomial  of 
degree  ♦  2-  in  j>  . 

P  given  by  (4.13)  has  generally  to  be  computed.  However, 


The  tensor 
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for  illustration,  it  will  be  evaluated  for  a  spherical  inclusion  in  an 
isotropic  matrix.  Thus,  L  fs  taken  to  have  components 


*  S  Cm 


*  ^  5* "  r 


(4.14) 


so  that  it  has  bulk  modulus  K-  and  shear  modulus  JJL  .  It  is  helpful  to 
use  the  symbolic  notation 


=  (iK>  2/*) 


(4.15) 


introduced  by  Hill  (1965a),  in  which  the  product  of  two  isotropic  ten- 

5  A  =  and  is  Rivt*n  by 


sors 


a  |  =  (eO(jO,(v>HvO)- 


(4.16) 


It  is  easily  deduced  from  (4.14)  that 


A< 

r-  p 

and  the  components  of  |  follow  correspondingly.  The  tensor  J  is 
isotropic  and  its  'moduli'  tC_ 


(4.17) 


>  ’A 


mav  be  obtained  from  the  scalars 


?...  =  Ifc. 


«kk 


,  P ...  , 

V'J 


IK,  *  (OAy 


(4.18) 


for  each  of  which  the  integrand  of  (4.10)  reduces  to  a  constant.  Com¬ 


pleting  the  details  gives 
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(  f — — 

I  3AC  + 


4/1  SyU  (3/C+  4-/0 


SfiL  ). 

+  4J0S 


(4.19) 


No  other  special  case  is  particularly  simple.  However,  some  reduc¬ 
tion  is  possible  when  the  inclusion  is  long  and  thin,  so  that  it  approx¬ 
imates  to  a  cylinder  with  elliptic  cross-section.  If  its  principal  axes 
are  aligned  with  coordinate  axes  and  the  3-axis  is  the  one  of  greatest 
length,  the  matrix  A  may  be  taken  diagonal,  with  entries 
where  £  is  small.  The  integral  in  (4.13)  may  now  be  transformed  by 
projection  onto  the  circular  cylinder  of  unit  cross-section  defined  by 

.  The  element  of  area  on  this  cylinder  is 
i4*\,  ,  representing  arc  length  around  a  cross-section,  and  the  ele¬ 
ment  of  solid  angle  4S  is  given  by 


iS  did j, 


iVJ/* 


(4.20) 


The  function 


no 


is  unchanged  by  the  projection  because  it  is  homo¬ 


geneous  of  degree  zero.  The  variable  is  now  replaced  by  K‘U* 


and  (4.13)  takes  the  form 


P  =  a,a 


.  1 4.  U;  r(i)[*;iyx>zr/lo>  *T):’A 

iK \x-  i  J  . 


K'\* ' 


(4.21) 


The  integral  with  respec 


tto  i: 


may  be  estimated  asymptotically  by 


setting  £  s  D  directly;  performing  the  resulting  elementary  integration 
with  respect  to  then  gives 
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(4.22) 


p~  v?£  * 

i,r  V*c-<  ~ 

Equation  (A. 22)  offers  a  solution  to  the  cylindrical  inclusion  problem 
that  is  an  alternative  to  the  usual  one  involving  complex  variables. 

Reverting  now  to  composites,  suppose  that  a  composite  comprises  a 
matrix  with  moduli  L  »  which  contains  rt  different  types  of  ellipsoidal 

M  * 

inclusion,  with  moduli  L  ,  O  ,  distributed  at  volume  con¬ 
centrations  C.  ,  where  each  is  small.  The  overall  moduli  of  the 

c  r 

composite  may  be  estimated  by  imposing  a  uniform  strain  fi.  at  infinity. 

A* 

The  mean  stress  O"  in  the  composite  is  given  by  the  formula 


n 

^  -  L  ^  ^  Crlr  , 

AM  ^  ^  .1 


(A. 23) 


where  represents  the  mean  polarization  in  the  r*th  phase.  To  lowest 

order  in  C ^  ,  the  inclusions  do  not  interact  and  may  be  estimated 

as  the  solution  of  equation  (A. 9),  with  ■£  =  ,  SUL  -  [_  and 

P  •*.  P  ,  the  value  appropriate  to  an  inclusion  of  rth  type.  The  low- 

Aw  ••  f 

concentration  estimate  of  the  tensor  of  overall  moduli  then  follows > 
using  (2.1),  as 


l  ~  L  +  Z- 

-  -•  r.i 


a^-L.y*  pj". 


(A.2A) 


This  result  is  (asymptotically)  exact.  It  is  affected  by  the  symmetry  of 
the  individual  inclusions  but  not  by  their  relative  positions  because 


5. 


THE  SELF-CONSISTENT  METHOD 


Consider  now  a  composite  comprising  a  matrix  containing  n  dif¬ 
ferent  types  of  inclusions,  not  necessarily  ellipsoidal,  with  tensors  of 
moduli  L  and  volume  concentrations  ,  f*  s  I  X'"  rv  ,  not  neces- 
sarily  small.  A  comparison  material  different  from  the  matrix  will  be 
used,  so  the  matrix  is  now  regarded  as  phase  tl  ♦  I  ,  with  moduli 

A4I 

and  volume  concentration  C_  •  Let  the  average  values  of  the  stress 

n4 1 

and  strain  over  the  ^th  phase  be  denoted  <ST  ,  .  Then 

r  —•  f 


■r  *  kr  tr  , 


(5.1) 


and 


r=.i 


cr*.r 


(5.2) 


*♦1  . 

<T  = 


A*l 

z. 


Eliminating  C  t  between  (5.2)  and  (5.3K  gives 

*1+1  *+  “I  z 

n 

W  ■  L  C-  +  ^  r  kr”  kfl+i  /  ~  r 

•*-  A+l  ***  rai  r  r 


(5.3) 


(5.4) 


so  that,  if  £.  is  prescribed  and 


t  e  A  ft. 


€. 

r  /v%>  < 


the  tensor  of  overall  moduli  L.  is  given  by 


(5.5) 
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k  -L.,* 

A.. 


(5.6) 


The  problem  now  is  to  estimate  the  tensors  At  .  An  approximate  method 
for  performing  this  estimation  is  to  embed  a  single  inclusion  of  rth  type 
in  an  infinite  matrix  composed  of  the  overall  material.  If  the  inclu- 
is  an  ellipsoid,  A  follows  from  the  results  of  the  preceding 


sion 


L  =  L  and  SL-L  -L_. 

"V*  O 

found  from  equation  (4.9)  with  and  then  equation  (3.9)  gives 


section,  with  L  a  L  and  iL.®  L  ~  L.  .  The  polarization  *£  is 

ZLo  SZ  —  -v P 


(5.7) 


If  the  Pth  inclusion  is  not  an  ellipsoid,  the  solution  of  the  inclusion 
problem  is  less  straightforward.  It  is  convenient  for  the  discussion  to 
follow,  however,  to  treat  (5.7)  as  defining  P  in  such  cases.  Sub- 
stituting  (5.7)  into  (5.6)  now  gives  an  algebraic  equation  for  L  .  since 
it  appears  on  the  right  side  both  directly  and  implicitly  through  the 


tensors 


—  r 


The  estimate  for  L  that  results  from  solving  this  equa¬ 


tion  is  usually  termed  the  self-consistent  estimate.  The  prescription 
was  developed  in  essentially  this  form  by  Hill  (1965c)  and  Budiansky 
(1965).  It  allows  for  inclusion  shape  and  also  approximately  for  inter¬ 
actions  because  the  inclusion  is  embedded  in  the  'overall'  material.  It 
also  reduces  correctly  to  (4.24)  in  the  limit  of  low  concentrations.  A 
possible  objection  is  that  the  inclusions  actually  interface  with  'over¬ 
all'  material  rather  than  with  the  matrix.  One  method  of  overcoming  this 
is  to  surround  the  inclusion  by  a  'shell'  of  matrix  material  and  then  to 
embed  this  new  compound  inclusion  in  an  Infinite  matrix  of  'overall' 
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material.  For  spherical  inclusions,  it  is  reasonable  to  take  the  shell 
as  spherical  and  to  choose  its  radius  so  that  the  inclusion  concerned 
occupies  exactly  its  share  of  the  matrix.  Kerner  (1956)  solved  this 
problem  for  isotropic  phases,  exactly  for  an  imposed  dilatation  and 
approximately  for  an  imposed  shear.  The  exact  solution  for  an  imposed 
shear  has  been  provided  by  Christensen  and  Lo  (1979). 

Suppose  now  that  the  composite  has  no  clearly  defined  matrix  phase 
or,  equivalently,  that  ®  .  The  prescription  (5.6)  fails  in 

this  case  but,  if  the  shapes  of  the  phases  can  still  be  distinguished, 
it  remains  possible  to  estimate  the  tensors  $  by  solving  inclusion 
problems.  Three  possible  prescriptions  for  estimating  L.  are  now  avail- 
able.  Equation  (5.2)  requires  that 


T  \  1"  I 


•  I  . 


(5.8) 


equation  (5 . 3) 2  requires 


(5.9) 


and  finally 


,  ©•*  L  £  *21cf  T 

mm  ^  «•*  mmW 


.implies 


=  0 


(5.10) 


Equations  (5.8)  to  (5.10)  are  not  all  equivalent,  for  (5.8)  can  be  ex¬ 
pressed  in  the  form 

-I 


ui.  nt,-D“ar  -  * 


(5.11) 
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and  equation  (5.10)  is  equivalent  to 


i  •- 1, 1 1  •  t  (L.-  Df  If !,  c<  fi  *  £.  ( t,-  DI'T'- 


(5.12) 


It  is  clear  from  these  relations  that  equations  (5.8)  to  (5.10)  become 
equivalent  when  f  s  P  for  all  r  ,  that  is,  when  all  phases  can  be 
modelled  by  similar  ellipsoids.  The  prescriptions  (5.8)  to  (5.10)  were 
proposed  variously  by  Hershey  (1954),  Kroner  (1958)  and  Eshelby  (1961) 
for  estimating  the  overall  moduli  of  a  polycrystal  by  treating  each 
grain  as  a  sphere.  Their  equivalance  was  proved,  essentially  as  above, 
by  Hill  (1965a).  It  is  true  also  that,  when  *P  s  F  for  all  f*  , 

**  f  •%« 

equation  (5.6)  can  be  manipulated  into  any  of  the  forms  (5.8)  to  (5.10), 
with  n  replaced  by  Of*  I  .  This  was  noted  for  spheres  in  an  isotropic 
matrix  by  Budiansky  (1965)  and  generally  by  Hill  (1965c). 

The  great  advantage  of  the  self-consistent  method  is  its  simplicity. 
Furthermore,  it  will  be  seen  when  examples  are  discussed  later  that  it 
often  yields  results  of  acceptable  accuracy.  There  are,  however,  some 
situations  with  which  it  cannot  be  expected  to  cope.  If,  for  example,  an 
isotropic  matrix  contained  spherical  isotropic  inclusions,  arranged  on  a 
lattice,  the  simple  self-consistent  method  would  estimate  the  overall 
moduli  as  isotropic  even  though  they  should  plainly  reflect  the  symmetry 
of  the  lattice.  Although  exact  methods  are  available  for  periodic  com¬ 
posites,  real  composites  are  usually  not  perfectly  periodic,  though  they 
may  well  display  some  directional  bias  'on  average';  the  inadequacy  of 


! — Jl 
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the  self-consistent  method  would  remain  in  this  case  also.  The  section 
that  follows  presents  a  method  that  makes  systematic  allowance  for  such 
phase  geometry. 
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6.  BOUNDS  FOR  OVERALL  MODULI 

If  every  detail  of  the  geometry  of  a  composite  were  known,  overall 
moduli  could,  in  principle,  be  calculated  exactly  by  solving  the  rele¬ 
vant  boundary  value  problem.  In  practice,  however,  except  in  cases  such 
as  those  displaying  periodicity,  a  complete  solution  could  not  even  be 
computed  because  of  the  detail  that  it  would  contain.  Also,  in  practice, 
only  certain  'average'  features  of  the  microstructure  will  be  known.  If, 
for  example,  only  the  moduli  and  mean  concentrations  of  the  phases  are 
known,  the  only  completely  accurate  statement  that  can  be  made  is  em¬ 
bodied  in  the  inequalities  (2.15);  that  is,  that  the  Reuss  and  Voigt 
estimates  bound  the  overall  moduli.  If  more  information  is  available, 
however,  it  should  be  possible  to  provide  better  bounds.  This  inform¬ 
ation  is  likely  to  be  statistical  in  character  and  it  is  necessary  now 
to  make  a  slight  digression  to  discuss  random  media  and  introduce  there¬ 
by  the  quantities  that  will  be  employed. 


6.1  Random  Media 

A  random  medium  is  one  of  a  family,  any  member  of  which  may  be 
characterized  by  a  label  &  that  belongs  to  a  sample  space  $  .  The 


value  of  0 C  is  taken  as  defining  the  medium  completely.  It  has  asso¬ 
ciated  with  it  a  probability  density  fW  which  is  defined  over  ^  . 
The  sample  space  ^  could  be  of  low  dimension,  for  example  if  the  ran¬ 


dom  medium  were  known  to  be  perfectly  periodic,  except  for  uncertainty 


over  the  position  of  any  one  representative  cell,  or  it  could  be  of  in¬ 
finite  dimension,  with  representing  the  values  of  f  at  each  indivi¬ 


dual  point  of  the  medium.  For  an  0 -phase  material,  it  is  convenient  to 
introduce  the  indicator  function  l(*)  ,  that  takes  the  value  1  if 
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3t  lies  in  phase  f*  and  zero  otherwise;  it  depends,  of  course,  on  c<  . 

The  mean  value,  or  ensemble  average,  of  defines  the  probability 

of  finding  phase  C  at  X  .  Thus, 

JS 

Likewise,  the  probability  P  (x  X f)  of  finding  simultaneously  phase 

r  5  v  / 

0  at  X  and  phase  J  at  X/  is 

f,,  (i ,1')*  (fjs)l(z'))  -  I  <6-: 

Probabilities  involving  more  points  are  defined  similarly.  If,  now, 
phase  r  has  moduli  L  .  ,  0  el,  2"’  ft  ,  the  value  L.  ( of  the 

<  /  MV  / 

tensor  of  moduli  at  X  can  be  written 


(6.2) 


Uj)  *  X  LJJ*-) 


(6.3) 


and  this  has,  using  (6.1),  mean  value 


<k(*)>  -  ZkrU*) 


(6.4) 


Similarly, 


r*  «■» 


<L(*)©Lt*o>  - x i  Lr©Lf 


(6.5) 


The  individual  components  in  (6.5)  represent  two-point  correlation  functions. 
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Correlation  functions  for  more  points  follow  a  similar  pattern;  clearly, 
they  are  expressible  in  terms  of  i Lf*.  X^")  and  so  on,  for 

any  physical  quantities  whose  values  are  given  for  each  phase. 

The  foregoing  description  is  useful  for  developing  statements  ’on 
average',  taken  over  many  specimens  produced,  perhaps,  by  some  manu¬ 
facturing  process  with  which  the  probability  is  associated.  A 

statistically  uniform  material  is  one  for  which  the  probabilities  w 
etc.  are  insensitive  to  translations;  thus  Us)  reduces 
to  a  constant,  becomes  a  function  of  (&-£')  only,  and 

so  on.  For  such  a  medium,  it  is  usual  to  make  an  assumption  of  ergodic 
type,  that  local  configurations  occur  over  any  one  specimen  with  the 
frequency  with  which  they  occur  over  a  single  neighbourhood  in  an  ensemble 


of  specimens.  In  this  case,  the  probability 


ility  P. 


is  identified  with  the 


volume  average  of  the  function  Us)  and  so  represents  the  volume 
concentration  cr  of  phase  r  ,  while  P,,  Ci,  £')  becomes  the  average 
of  k(2L+£*)fs  as  ranges  over  a  large  volume. 

Strictly,  these  notions  require  the  medium  to  be  infinite  but  they  pro¬ 
vide  useful  empirical  definitions  in  practice.  Clearly,  it  will  be 
possible  to  determine  multipoint  probabilities  only  up  to  some  relatively 
low  order.  Correspondingly,  Kroner  (1977)  proposed  classifying  materials 
as  uniform  ’of  grade  It  ’  if  probabilities  involving  up  to  It  points 
were  known  to  be  translation-invariant.  He  classified  other  properties 
(such  as  isotropy)  similarly. 

For  future  use,  the  two-point  probabilities  P  [  X  jU) 
discussed  for  a  two-phase  material.  They  satisfy  the  relations 


are  now 
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PjUQ  *  Us)  *  ■» 

?n fe >&')*■  &')  =  1,U)  ‘  s)  +  i), 

•li&'.s) ♦  (Us' s)  =  Us.)  =  £,<£,*')  *  ?n(*,*!)  <6-6> 


always.  Therefore,  when  the  medium  is  statistically  uniform  (up  to 
grade  two) , 


12,  1  f„=  c,c1(i-«))  Plt  =  <Uc,<U, 


(6.7) 


where  K  is  some  function  that  depends  upon  fx  — fcOonly  which,  further- 

-v  / 

more,  is  an  even  function  of  its  argument.  If,  in  addition,  the  medium 
has  no  long-range  order,  P  'v  C,  at  large  separations  of  X  and 
X7  and  tends  to  zero  as  jx-X^J  tends  to  infinity.  If  the  dis¬ 
tribution  of  phases  is  isotropic  (up  to  grade  Z  ),  then  "£  is  a  function 
of  |  X  —  x'  |  only;  it  should  be  noted  that  this  statement  makes  no  refer¬ 
ence  to  the  elastic  symmetries  of  the  phases  themselves  so  that  the  com¬ 
posite  could  still  be  elastically  anisotropic. 

Suppose  now  that  the  composite  comprises  a  matrix  in  which  are 
embedded  inclusions  of  known  shapes  and  orientations.  All  mean  values 
for  this  composite  follow  from  knowledge  of  the  probability  of  finding 
any  specified  number  of  inclusions  in  any  specified  region.  This  will 
not  be  discussed  in  detail  but,  for  purposes  of  later  application,  if 
the  matrix  contains  just  a  single  population  of  inclusions,  centred  at 
X.  As  I  2  ••••  ,  I  will  be  defined  as  the  probability  density  for 

>  >  A 
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finding  an  inclusion  centred  at  X.  ,  P  will  denote  the  joint 

—A  *0 

probability  density  for  finding  different  inclusions  centred  at 
X  ,  X  and  so  on.  Statistical  uniformity,  isotropy,  etc.  for  such  a 

~A  op... 

composite  is  expressed  in  terms  of  etc.  in  a  similar  manner 

to  that  described  above  for  P^  ,  etc. 

6.2  Hashin-Shtrikman  Bounds 

In  most  applications,  it  is  unlikely  that  information  of  grade  higher 
than  two  (in  Kroner's  terminology)  will  be  available.  Therefore,  although 
formal  extensions  are  possible  (and  will  be  outlined  later),  attention 
will  be  focussed  now  on  the  best  bounds  that  can  be  obtained  using  two- 
point  correlations  only.  Inequalities  such  as  (3.13)  which  are  derived 
from  the  classical  energy  principles  are  excluded  immediately  for,  written 
explicitly,  the  last  term  in  (3.15)  is 

V-  J<l*  f  ■» s'  IY« .  &')X*(a')*k<&) V 


which  involves  points  taken  three  at  a  time.  The  Hashin-Shtrikman  ine¬ 
qualities  (3.20),  (3.26)  were  constructed  to  contain  no  such  term  and 
bounds  can  be  obtained  from  them  by  considering  trial  fields  of  the  form 


(6.8) 


r-( 


for  an  fl  -phase  material , where  the  X  are  functions  of  x  alone,  inde- 
pendent  of  the  configuration:  any  relaxation  of  this  restriction  would 
introduce  terms  into  fit*)  that  would  involve  more  than  two  points, 
so  that  (6.3)  is  the  most  general  trial  field  that  can  be  allowed. 
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Substituting  (6.8)  into  (3.17)  now  gives 

f(l*)  -  V'ji  j's  .)'!,(*)  -  2  £.(*)] 

+L±Uih  U)L  ($WZU,  s')  l&)hU')j 

r *1  *ei  ;  </  J  (6.9) 


This  depends  upon  the  sample  ,  as  it  must.  Further  simplification  is 
possible  subject  to  the  assumption  that  the  functions  J"  vary  rapidly 
over  regions  that  are  still  sufficiently  small  relative  to  the  macro¬ 
scale  for  the  functions  %(*),[(*,  s')  to  be  treated  as  effectively 
constant.  It  is  credible  (though  not  rigorously  proven)  in  this  case 
that  the  value  of  HtV  is  independent  of  o<  and  so  identical  with 
its  mean  value,  or  ensemble  average.  Then,  (6.9)  reduces  to 

Ht*)‘  jdsPJi)Zr(i)[(LrLoyifCi)-iiJi)J 

I  f‘l  * 

*ti~  Us  IrU)!**'  f(Z,  i'jTsU')  %Ji,  &') 

r -|  J*|  J  J  ** 


(6.10) 


This  is  extremized  when 


*  £ .(*)  ;  fl,*-  n, 


(6.11) 
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to  give  the  extreme  value  —  (<1>,  &.'),<!) 
average  of  ^ 


denoting  the  ensemble 


It  may  be  noted  that,  if  the  integral  equation  (3.16)  is  averaged 
conditionally  upon  phase  C  being  present  at  X  ,  it  gives 


=  ^  O) 


(6.12) 


where  <1  (*)>,  denotes  the  expectation  value  of  Z(  5)  ,  conditional 

upon  phase  being  present  at  fC  //  denotes  the  expecta- 

***  '/rj 

tion  value  of  T  (x')t  conditional  upon  phases  r  and  S  being  present 
at  x  .  respectively  and  p*(  5',*)  is  the  corresponding  con¬ 
ditional  probability  for  finding  phase  S  at  Xf  ,  so  that 


L *  L  U',a)K(z) 


(6.13) 


Equation  (6.12)  reduces  to  (6.11)  if 


and,  in  addition 


<i>; 


<!(*')>[  =  <£(!')>!■ 


is  identified  with  t 


(6.14) 


Equation  (6.14)  is  not,  of  course,  strictly  correct,  unless  the  composite 
has  some  special  structure  such  as  periodicity.  It  is  an  example  of 
the  'quasicrystalline  approximation',  introduced  as  an  ad  hoc  approxi¬ 


mation  by  Lax  (1952);  its  appearance  from  the  variational  principle  at 
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least  partially  explains  Its  success. 

The  discussion  so  far  has  been  rather  general;  it  is  desirable  now 
to  specialize  to  the  case  of  a  statistically  uniform  medium  and  to  the 
boundary  condition  (2.2),  so  that  e.#  s  to  obtain  bounds  upon  over¬ 
all  moduli.  The  body  has  already  been  assumed  to  be  large  relative  to 
the  microstructure  and  it  might  be  expected  that  ?(*,  t')  could  be  re¬ 
placed  by  its  infinite-body  form  (4.5).  This  is  not  quite  correct  be¬ 
cause  this  is  of  order  I  i  |  at  large  I  X  J  so  that  the  integral  in 
(6.11)  would  converge  only  conditionally.  It  may  be  noted,  however,  that 
the  mean  value  of  the  exact  P*  over  \/  is  zero,  since  it  is  associated 
with  zero  boundary  displacement.  Consequently,  for  any  field  *?*  that 
oscillates  rapidly  about  a  mean  TT  ,  /  “T  ~  I  [  'C—  T / exactly  and,  in 
this  latter  form,  a  convergent  integral  is  obtained  if  |  is  now  re- 
placed  by  {  .  This  does  not  prove  that  this  prescription  is  correct, 

but  further  arguments  for  its  plausibility  have  been  advanced  by  Korringa 
(1973),  Willis  and  Acton  (1976)  and  Willis  (1977),  which  are  also  developed 
slightly  differently  in  the  review  of  Willis  (1980d).  These  essentially 
hinge  upon  representing  P*  in  terms  of  the  infinite-body  operator  as 
above,  plus  an  additional  term  whose  influence  will  be  felt  only  in  a 
layer  close  to  the  boundary  w .  In  this  approximation,  the  operator 
P*  is  insensitive  to  translations  and  the  functions  T  reduce  to  con- 

/w  . 

stants  Y  that  satisfy  the  algebraic  equations 


(6.15) 
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where 


V  ( i') 


represents  the  probability  of  finding  phase  5  at 


jc'  ,  given  phase  T  at  the  origin.  The  concentration  C  is  its 

**  9 

value.  Having  solved  the  system  (6.15),  to  give 


mean 


X  S  e 

r  -•«  r  1 


say,  the  extreme  value  of 


Jfa*) 


is  expressible  in  the  form 


(6.16) 


(6.17) 


where  the  estimate 


L  =  L  i  cn 

AM  Aa  *'  ^  | 


(6.18) 


for  the  overall  moduli  bounds  L.  from  below  if  L  -  L  is  positive 

—  r  —  o 


L.  from  below  if  L  —  L. 

— r  — 

definite  for  each  f*  ,  and  from  above  if  L^.”  L.  #  is  negative  definite 
for  each  C  . 

The  solution  (6.18)  takes  a  particularly  simple  form  if  all  of  the 
functions  V  U')  are  isotropic  and  so  depend  upon  J  X'j  only. 

The  result  given  in  Section  4,  that  the  integral  of  I  over  any  sphere 
is  P  ,  independently  of  its  radius,  implies  that  the  integral  of  f 
over  any  shell  J**)  -  C  >  O  is  zero.  The  integrals  in  (6.15)  may 
therefore  be  replaced  by  integrals  over  a  small  sphere  centred  at  the 


origin,  within  which 


•sir 


differs  negligibly  from  its  value  at  the 


origin,  namely,  $  .  Evaluating  the  integral  then  gives 


**5 


(6.19) 


where  <*t  >*  2S  &• 


It  follows  now  that 
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(6.22) 
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Equations  (6.15)  now  take  the  form  of  (6.19),  with  n»2  and  P  re¬ 
placed  by  P  and  their  solution  yields  the  estimate  (6.21),  similarly 
modified,  for  L.  . 

Willis  (1977)  suggested  that  the  estimates  (6.21),  with  P  taking 
'ellipsoidal'  form,  might  be  appropriate  to  discuss  a  composite  containing 
aligned  ellipsoids  including,  as  limiting  cases,  thin  needles  and  flat 
cracks.  Justification  was  provided  by  Willis  (1978),  through  explicit 
study  of  a  composite  containing  inclusions,  whose  distribution  was  de¬ 
scribed  by  the  probability  densities  f,  ,  etc.  for  the  positions 

A  B|A 


of  their  centres.  If 


PlA 


was  'ellipsoidal',  (6.21)  was  reproduced. 


Willis  (1980a)  also  generated  a  result  rather  like  the  one  above  with  P 
replaced  by  P  ,  for  a  matrix  containing  a  single  population  of  aligned 
inclusions  of  any  shape,  distributed  arbitrarily.  Only  the  'ellipsoidal' 
example  has  been  worked  out  in  detail,  however. 

6.3  Bounds  for  Compliances 

<v 

Bounds  for  compliances  are  obtained  by  inversion  of  L  ,  to  produce 

/V  - ’  ** 

corresponding  estimates  M  for  M  .  These  can  also  be  obtained  directly, 
however,  by  considering  the  traction  boundary  value  problem,  with  boundary 
condition  (2.5).  The  bounds  follow  by  extremizing  the  functional  j,  (f) 
given  by  (3.30):  their  structure  is  exactly  like  that  for  L  ,  except 

that  ( L  —  L  )  is  replaced  by  (M  —  M)  and  f*  is  replaced  bv  •“  Ci 

— »*  o  —  I 

The  infinite-body  form  for  the  operator  ^  can  be  deduced  bv  noting 

I 

that,  when  V  is  large,  equations  (3.16)  and  (3.29)  must  be  consistent, 
with  T  i  L  ,  C  s  C.  and  —  «*"  —  [_  £.  f  V  The  result  is  that , 

**  -»p  *  "*»  n,|)  a,  '  —  ' 

for  large  V  ,  A.  is  defined  bv 

'*•  I 


&, :  7 


•v%.  ^ 


(6.23) 


*2* 
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Then,  for  a  composite  with  "ellipsoidal'  symmetry,  for  example. 


■  |Z  ■ Cl  >z  «*  U  «.)£  fa, , 

1  (6.24) 


where 


L  P  L  -  L 


«w  o  *•  a 


x:  o 


(6.25) 


/s/  /V 

It  can  be  verified  directly,  using  (3.24),  that  !_  and  1*) ,  as  given 

•O' 

by  (6.21)  and  (6.24),  are  inverses. 

6.4  Bounds  of  Higher  Order 

The  bounds  that  are  given  by  (6.18)  are  not  the  best  that  are  ob¬ 
tainable  from  the  trial  field  (6.8)  for  *  since  they  were  derived 
from  the  Hashin-Shtrikman  principle  rather  than  the  stronger  classical 
principles.  The  penalty  for  using  the  latter  is  that  three-point  proba¬ 
bilities  will  also  be  involved.  Still .assuming  that  these  are  known, 
(6.8)  is  now  substituted  into  the  minimum  energy  principle  (3.15),  to 
give 

V~'jdx  I  t,(x-)Z  (L'-LjlUjeJx) 


(6.26) 
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where 


r 

-^5 


*,(*.)  jdz'F  (*> 


(6.27) 


is)  -*  ^  r«'  xJirQI^/fe'T  (s,  i'lljx'JUij. 


((■a) 


Specializing  now  to  a  statistically  uniform  medium,  the  integrals  that 
appear  in  (6.27)  and  (6.28)  may  be  evaluated  as  in  the  preceding  section: 
for  (6.27),  (*0  is  replaced  by  fj( X'J  “  C(  and  f~*  is  replaced 

by  P*  ,  and  similar  replacements  are  made  in  (6.28).  The  expectation 
values  of  the  operators  A  »  8  then  become  translation-invariant 

-v  fj  — 

and,  if  &■  f  X]  is  taken  equal  to  &  ,  the  functions  T  may  be 

taken  constant.  If  the  right  side  of  (6.26)  is  now  replaced  by  its  ex¬ 
pectation  value,  it  follows  by  extremizing  that 


I  Li  =r  sUs-Xf  ZM-XL-Ut, 

where  L.  is  the  Voigt  estimate  and  T  /*  =•  f.  2  **•  fl 

a.  y  •**"  I  ' 

satisfy  the  equations 

Z  (<!,>  ■*<!„}) 5.  -  ? <L>(L-L,)l , 


(6.29) 


(6.10) 
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with 

<AM>  '  J-5'f  Cj)  <6-31> 


and 


Following  Beran  and  Molyneux  (1966)  and  Dederichs  and  Zeller  (1973), 
Kroner  (1977)  generated  an  'optimized'  third-order  bound,  essentially  by 
substituting  into  (3.15)  the  trial  field 


X*  =  fi-  t,  ,  (6. 

in  which  the  constant  tensor  was  arbitrary.  This  possibility  is 

included  within  (6.8)  and  the  bound  (6.29)  is  correspondingly  tighter. 

The  work  required  to  produce  the  Beran-Molyneux  bound  is  only  marginally 
less  than  is  needed  to  produce  (6.29),  since  it  requires  evaluation  of 

<ftr*>and 

Kroner  (1977)  also  produced  a  tighter  bound  for  a  special  type  of 
material  that  he  termed  'disordered'.  Its  distinctive  feature  (which  he 
describes  in  terms  of  the  moduli  and  topology  being  uncorrelated)  is  that 
it  permits  the  replacement  of  the  awkward  function  ft£)  in  (6 . 3j)  and 
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(6.32)  by  £«S>  ,  where  the  constant  tensor  P  is  defined  by  (4.10). 
For  (6.31),  this  is  legitimate  so  long  as  F  is  isotropic,  but  for 


(6.32)  it  requires  more  justification:  Kroner  (1977)  considers  the  re¬ 
placement  justified  so  long  as  |_  is  identified  with  1.  .  ,  when 


—  • 


.  Without  arguing  for  or  against  this  prescription,  it  is 
simply  remarked  that  its  adoption  yields  simple  forms  for 
for  which  (6.29)  and  (6.30)  lead  to  the  upper  bound  (6.21),  previously 
derived  from  the  Hasliin-Shtrikman  principle,  except  that  now  is 

identified  with  L  .. 

—  V 

Analogous  lower  bounds  of  third  order  can  be  obtained  from  (3.22); 
no  explicit  development  will  be  given. 

Bounds  that  involve  still  more  statistical  information  can  be  found 
by  allowing  trial  fields  more  complex  than  (6.8).  For  example,  the  trial 
field 


l*(l)  *  U*)[lr  (i)  »  fjfr,  (*•&%&)**'], 


(6.34) 


in  which  the  functions  ~r  U)  and  frs  (*,£')  are  to  be  determined 
optimally,  introduces  correlations  involving  up  to  four  points  in  the 
Hashin-Shtrikman  principle,  and  up  to  five  points  in  the  classical  prin¬ 
ciples.  These  are  not  considered  further,  except  to  remark  that  Willis 
(1978)  employed  a  polarization  similar  to  (6.34)  for  a  matrix  containing 
a  single  set  of  inclusions  and  obtained,  by  considering  only  the  terms 
of  low  order,  a  variational  estimate  for  the  C*-  coefficient  in  the  low- 
concentration  expansion  for  L  ,  C.  being  the  volume  concentration  of 


the  inclusions. 
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6.5  Further  Comments  on  the  Self-Consistent  Method 
The  Hashin-Shtrikman  estimate  (6.18)  for  L  provides  a  lower  bound 
when  L  is  sufficiently  small  and  an  upper  bound  when  L  ^  is  suffi- 
ciently  large.  At  intermediate  values  of  ,  it  does  not  yield  a  bound, 

but  it  is  still  a  variational  estimate  and  it  might  be  postulated  that 


this  estimate  will  be  closest 


to  L 


when 


L  =  L  , 


for  then  the  mean 


polarization  will  be  zero.  An  estimate  of  self-consistent  type  is  there¬ 
fore  obtained  by  postulating  that  L.  is  identical  to  1_  when  /,  s  L  , 
That  is. 


•v  O 


L(L)  -  L 


(6.35) 


or,  equivalently. 


2.  ( \=0  =  L ) 


(6.36) 


A/ 

The  prescription  (6.35)  of  course  is  not  exact  because  L.  is  estimated 
from  a  piecewise-constant  polarization  field.  If  the  composite  has 

A/ 

isotropic  (or,  more  generally  'ellipsoidal')  phase  geometry,  L  is  given 
by  (6.21)  and  it  is  easily  seen  that  the  self-consistent  prescription 
(6.35)  yields  any  of  (5.8),  (5.9)  or  (5.10),  with  P  s  P  for  all  P  . 
The  interesting  point  is  that  the  prescription  (6.35)  allows  for  two- 
point  statistics  in  a  general  way  and  makes  no  reference  to  particular 
inclusion  shapes.  Thus,  it  provides  some  justification  for  the  use  of 
(5.8)  to  (5.10)  for  composites  such  as  isotropic  polycrvstals,  in  which 
the  grains  are  certainly  not  all  spheres,  and,  in  addition,  offers  an 
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extension  of  the  self-consistent  method  to  composites  in  which  the 
microstructure  displays  spatial  bias. 
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7  EXAMPLES 

The  application  of  the  preceding  theory  will  now  be  illustrated  for 
some  simple  examples. 

7.1  An  Isotropic  Mixture  of  Isotropic  Phases 

In  the  notation  of  Hill  (1965a),  introduced  at  equation  (4.15),  let 
a  composite  have  O  phases,  with  moduli 

distributed  statistically  uniformly  and  isotropically.  It  is  reasonable 
to  take  the  comparison  material  as  isotropic  so  that  L,  =  (3k,  z/<)  , 
say.  The  Hashin-Shtrikman  bounds  for  this  composite  then  follow  from 
(6.21),  with  the  tensor  ^  given  by  (4.19).  Explicitly, 

£  =  /  Z  v  l  Z  «,  fc.  , 

£  r  C3fcr*4u))  s 


/V 


Jic,  - f  Z  c,  _  {^1K^r),h  . 

t  *  t>J*f  (K  +Zp)  +/*  i  ?/))  *  (fs  (&+*/)  *  2f) 

(7.1) 


Equations  (7.1)  yield  upper  bounds  for  /c  whenever  /C  are 

chosen  larger  than  for  each  C  and  the  smallest  upper  bounds 

are  obtained  by  setting  /C.  C^r')  >  ~  (/*t)  ' 

Similarly,  the  greatest  lower  bounds  are  found  by  setting  /C«-  z  Mi* 

-  Mi*  fa)  .  These  bounds  were  given  by  Hashin  and  Shtrikman 
(1963),  with  the  implicit  restriction  that  had  both  to  be 

obtained  from  the  same  phase,  and  ’/V  similarly.  This  restric¬ 

tion  was  removed  by  Walpole  (1966a). 


Equations  (7.1)  also  provide  self-consistent  equations  for  the 
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determination  of  fC  ,  ,  if  K.  and  fc.  are  identified  with  /c.  and 

— 

and  are  identified  with  .  For  a  two-phase  material,  elimination 

of  K-  from  equations  (7.1)  yields  a  quartic  equation  for  .  This 
will  not  be  given  but  a  difficulty  with  the  self-consistent  method  will 
be  highlighted  by  considering  the  special  case  a ^tA.t  —  O  ,  so  that 
the  composite  represents  a  porous  medium  whose  matrix  phase  has  moduli 
fC L,  .  It  is  easy  to  show  in  this  case  that 


K  =  , 

3c,  kz  ♦ 

while  J*  satisfies  the  quadratic  equation 


(7.2) 


It/*  =  o. 

(7.3) 

If,  now,  ,  equations  (7.2)  and  (7.3)  have  solutions 

c,  O-c,) 

~  £  0  ~  ci)/*  i 

(3-S)  (7>4) 


The  corresponding  Hashin-Shtrikman  upper  bounds  are  found  from  (7.1) 
with  fC  ■  -Z : 


c  C*~  ^  fy*xO  c^) / •^^'0 


(7.M 


The  Hash  I  n-Slit  r  Ikm.in  lower  bounds ,  obtained  willi  /C  —  £P 


.uc  ,'i'in. 


K.|uat  Ions  (7. ''il  and  (7.S)  .iv;nr  at  low  values  ol  lull  I  lie  sell 

consistent  estimates  (7.'i)  reduce  I  o  :eio  al  a  pole  concent  I  al  1  on  e.|uul 
1 1'  1/7.  This  is  unreal  istie  and  demonsl  tales  I  lie  need  loi  caul  ion  in 
applvinp  the  se  1  i-cons  I  s  t  enl  method  In  extreme  cases.  Hill  (I 'hi  Sc)  and 
Hudlanskv  (1'lti‘i)  both  discussed  the  example  pteseiite.l  in  this  seel  ion 
with  0=1  ,  In  the  context  ol  a  matrix  cent  a i n i up  spherical  inclusions 

but,  with  l  he  interpretation  ol  Section  <>.‘>  .  auv  phase  p.eomet  i  v  is 
allowed,  si'  lonp  as  li  is  isotropic. 

Kor  a  two-phase  compos  1 1  e  tot  which  ate  both  .lei  ive.l  I  ioiii 

one  phase  and  ,  A  are  both  derived  Item  the  ol hei ,  the  bounds  oh 

tallied  from  (7.1)  apt  ee  with  the  estimates  ol  ovetall  moduli  piven  bv 
Kernel'  (ll,'i(i),  which  were  obtained  bv  sell  consistent  Iv  embed. line,  into  a 
matrix  ol  'overall'  material  a  composite  inclusion  conststuip  ol  a  sphere 
ol  one  material  encased  in  a  spherical  shell  ol  the  otbei,  the  l  ad  i  i  ho  hip. 
chosen  to  p.lvo  the  correct  volume  rat  ios  ol  each  phase.  l  itliei  bound 
can  be  reproduced,  bv  1  lit  erchanp,  I  Up.  the  two  phases.  ilit  islenson  and  I  o 
(ll)7'l)  noted  that  Kernel  's  solut  ton  lot  an  imposed  sheai  violated  a  con 
tilliiftv  condition,  t  lioup.h  his  solution  tor  imposed  dilation  Ol  compi  es 
Stoll  was  correct  .  Thev  provided  the  tine  solution  lot  the  she.it 
case  and  demonstrated  that  It  estimates  the  sheai  modulus  as  letup, 
between  the  Hash  I  n  Slit  r  i  km. in  bounds.  l-'oi  ext  tome  ditteiences  between 
phase  moduli,  it  mav  be  remarked  that  the  model  ol  Kei  net  and  ilii  Ist.'iiseu 


and  Lo  Ls  not  necessarily  more  reliable  than  tlio  simple  self-consistent 


model.  In  surrounding  an  inclusion  by  a  uniform  shell  of  matrix,  it  is 
likely  to  underestimate  interact  ions  between  inclusions  while,  to  judge 
from  (7. A),  the  simple  self-consistent  model  appears  to  overestimate 
them. 


7 . 2  1’ojvc  r vs  t  a  1  s 

A  polycrvstal  may  lie  regarded  as  a  limiting  ease  of  an  ft  -phase 
material  in  which  ft  — ^  a*  ,  each  phase  being  associated  with  a  different 
crystal  orientation.  Relative  to  crystallographic  axes,  any  one  crystal 

has  moduli  L  ,  sav ,  and  the  moduli  L  of  phase  C  are  obtained  from 

»t  r 

L  hv  a  rotation. 


c 

Consider  now  an  isotropic  polycrvstal,  in  which  each  crystal  orien¬ 
tation  is  equally  likely  and  all  of  t ho  functions  P  are  isotropic. 

list i mates  for  the  overall  moduli  L_  are  obtained  from  (h.Jl)  with 

■ 

taken  isotropic.  The  manipulations  are  easier  to  perform  if  (b.21)  is 
written  in  the  alternative  form 


where  <C  X  ^  =  Cf'  Xr 

I-  =  t(kr-La)'' <•£]*' 


Since  all  orientations  are  equal  lv  likely,  the  tensor  ^  is  i  so- 

t  t  op  ie ,  so  t  hat  <I>  =  .  'flu'  bulk  and  shear  'moduli' 

assoc i at  od  with  <I>  can  be  found  from  its  scalar  invariants,  as  given  hv 
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(4.18)  for  the  Isotropic  tensor  P  .  The  corresponding  scalars  calcu- 
lated  for  are  independent  of  rotation  and  so  may  be  calculated  for 

any  one  orientation,  say  that  for  which  L.^  =  Lc>  0nce  an<*  /*r 

have  been  found  in  this  way,  (7.16)  gives  the  bulk  and  shear  moduli  £ , 
J*.  in  the  form 

<  -  *  -*  , 


*  Mr  /( *  ~  , 


(7.8) 


where  L  (3/Cy2^tj.  If  Lt  is  cubic, 


so  that 


aw  ■  Sk  /  '  *“)  <  "Tx  -  V*  *■»  < 


(7.9) 


with  equations  for  other  components  obtained  by  cyclic  permutation  of  the 
suffixes,  it  follows  that 


tW)iiu  ’  **«  .  *  lK<  *  *  */•'  ■ 

Then,  if  the  notation  L  t  =  (iKCt  Z.Ja.1 , ^  introduced  by  Walpole 
(1966b)  is  employed,  products  and  inverses  can  be  worked  out  directly 


(7.10) 


T 

-  c 


T 


referred  to  the  cube  axes. 


and,  if 


represents 
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I.-  ( fc)  •  bfay [£-,  -a]"] 


Then,  from  (7.10)  and  (4.18), 


T  *  f  — 

l  S(K 


(7.12) 


Bounds  for  K.  follow  by  specifying  *’/  appropriately.  The  unique 

choice  for  fC  is  K.c  ,  which  gives  /Cj.  ®  0  and  /C  “  /Ct  ,  exactly,  , 

✓ 

on  the  other  hand,  may  take  either  of  the  values  A  *  A  •  the  greater 
giving  an  upper  bound  and  the  lower  a  lower  bound  for  /*  .  When s yUc  , 


/•A 


where 


'>A‘ J6 


and,  when 


/  ^ 


(7.13) 


,  where 


A  ■  x  •  *  for 


(7.14) 


In  these  equations,  ^(4  and  are  obtained  from  P  ,  evaluated  with 


and 


X 
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and  »  respectively.  The  bounds  (7.13),  (7.14)  were  first 

given  by  Hashin  and  Shtrikman  (1962b). 

A  self-consistent  estimate  may  be  obtained  from  equation  (6.23), 
which  implies  that  =  0  .  The  first  of  these  gives  AC  *  , 

exactly,  while  the  second  can  be  manipulated  to  give 


VO/1"  X  °  (7.15 

This  result  was  first  given  by  Hershey  (1954)  in  the  form  of  a  quartic 
equation  with  factors  (7.15)  and  (  *  1kc  ).  He  assumed  spherical 

grains  and  applied  the  prescription  (5.9). 

Anisotropic  polycrystals  require  a  computation.  Kneer  (1965)  studied 
textured  polycrystals,  in  which  grains  were  modelled  as  spheres  but  not 
all  crystal  orientations  were  equally  likely.  This  required  the  evalua¬ 
tion  of  P  for  an  anisotropic  comparison  material,  using  (4.10).  Ex- 
amples  for  which  the  actual  phase  geometry  is  anisotropic  do  not  yet 
appear  to  have  been  studied,  though  the  development  of  Section  6  now 
makes  this  possible. 

7.3  A  Composite  Containing  Long  Aligned  Fibres 

The  prototype  of  the  composite  now  to  be  considered  consists  of  a 
matrix  reinforced  by  long  parallel  fibres,  distributed  randomlv  on  any 
cross-section.  Such  a  composite  is  transversely  isotropic  and  was  first 
studied  systematically  by  Hill  (1964).  For  present  purposes,  there  is 
nothing  to  be  gained  from  immediate  specialization  to  this  extent  and  so  an 
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fl  -phase  composite  will  be  considered,  with  each  phase  transversely  iso¬ 
tropic  with  symmetry  axis  parallel  to  the  3-axis,  and  with  two-point 
probabilities  £  X  -  X'  )  that  depend  only  upon  C ~  )  , 

Greek  suffixes  taking  the  values  1,  2  only.  The  composite  may  thus  be 
treated  as  two-dimensional,  or  else  may  be  regarded  as  the  limit  of  a 
composite  with  'spheroidal'  symmetry.  In  either  case,  estimates  of  over¬ 
all  moduli  are  obtained  from  (6.21),  with  the  tensor  P  taking  the  form 
(4.16)  with  fltj  s.  . 

The  tensor  of  moduli  |_  is  most  usefully  expressed,  in  the  notation 
of  Hill  (1964,  1965b)  and  Walpole  (1969),  as 


L  -  2m,  2f) 


(7.16) 


when  <f~  -  l_  <£  i 


implies 


Si  “  ^ 

C*Ti  ~  *ii)  *  znfa r*it)  i 


%  "  2f  >  *31  =  Zf  €il 


(7.17) 


In  this  notation,  the  product  LI  of  two  such  tensors  is  given  bv 

I  I. 


57 


L,L=  (*> <7- 


(»/iA  , -t/2A,-  t'/iA,  k/a  ,  \/l~  ,  I /if), 


(7.19) 


where 


A  *  kft- W. 


(7.20) 


When  L  is  a  tensor  of  elastic  moduli,  it  has  the  symmetry  4  -  but 
products  do  not  share  this  property.  The  moduli  It  ,  m  are  the  plane- 
strain  bulk  and  shear  moduli.  Other  useful  quantities  are  the  plane- 
strain  Poisson's  ratio  >)  and  the  longitudinal  Young's  modulus  E  . 
These  are 


v  *  l/zk  ,  £  =  « -  r/* 


(7.21) 


Positive-definiteness  of 


positive. 


L  requires  that  all  of  k  ,  m  ,  E  and  j)  a 


Taking  the  comparison  material  as  transversely  isotropic  now,  witli 


components  as  in  (7.16)  with  , 

(4.16)  that 


it  is  fairly  easy  derived  from 


p  s  (\/z(k*M)t  0,  0,  +  5 


(7.22) 
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this  was  given  by  Walpole  (1969)  and  Laws  and  McLaughlin  (1978) .  All 
that  is  required  now  is  routine  substitution  into  (6.21).  Walpole  (1969) 
gives  the  results  in  the  form 

£  „  [£c./(fcf4*)J''-«  ,  t  = 

%  =  ^  X  C,  /fi*,  4  2«)]j  - 

f  *  '-f. 

(7.2 

from  which  both  bounds  and  self-consistent  estimates  follow.  Equivalent 
formulae,  together  with  useful  variants  involving  £■  and  V  ,  were  given 
by  Hill  (1964,  1965b)  for  a  two-phase  composite,  except  that  he  obtained 
no  bound  for  rw  .  llashin  (1979)  recently  reconsidered  the  two-phase 
problem  and  derived  estimates  of  moduli  from  a  'composite  cylinder' 
model  similar  to  the  spherical  model  of  Kerner  (1956).  It  generated  as 
estimates  the  bounds  that  are  obtained  from  (7.23). 

7.4  Composite  Containing  Short  Fibres 

Fibre-reinforced  composites  frequently  contain  fibres  whose  length 
to  thickness  ratio  is  of  the  order  of  10  to  100.  The  analysis  of  the  pre¬ 
ceding  section  therefore  applies  only  approximately  and  some  assessment 
of  its  range  of  validity  is  called  for.  Such  an  assessment  can  be  made 
if  the  fibres  are  modelled  as  aligned  spheroids.  This  allows  a  fairlv 
straightforward  self-consistent  analysis  and,  if  the  distribution  of  the 


1 


I 


i 


4 


I 


I 
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fibres  is  taken  as  ’spheroidal’,  bounds  can  also  be  constructed  using 
(6.21).  The  self-consistent  approach  requires  the  evaluation  of  f*  for 
a  transversely  isotropic  matrix.  Although  the  integrals  can  be  worked 
out  analytically,  their  form  is  still  sufficiently  complicated  to  require 
numerical  evaluation  by  computer  except,  of  course,  that  for  long  fibres, 
P  approaches  the  value  (7.22).  Laws  and  McLaughin  (1979)  have  carried 
out  the  required  computations  for  a  glass  reinforced  polyester  resin. 
Their  results  will  be  discussed  later.  A  complete  set  of  bounds  has  not 
been  found  but,  in  the  course  of  a  study  of  waves  in  composites,  Willis 
(1980b, c)  obtained  as  a  by-product  the  estimate  (6.21)  for  'spheroidal' 

P  ,  when  both  phases  were  isotropic  and  the  comparison  material  was 
identified  with  one  of  them.  Again,  general  results  could  be  obtained 
only  with  a  computation  but  som<j  general  feeling  for  the  influence  of 
fibre  length  can  be  gained  by  considering  the  asymptotic  form  taken  by  P, 
for  a  spheroid  whose  thickness  to  length  ratio  is  £  ,  where  £  «.  I  . 
Willis  (1980b)  gave  the  formula 


0^4/-) 


.  o,  0,  JK 


JtZtlt-  ,  —  ^ 


+  £l  U  (z/iQ_  ( -(} 

(-(**-/)>  }«■/>,  SKy*,  SKilfjl). 

’  V  '  (7.24) 


Russel  (1973)  did  no 


t  give  P 


but  derived  an  equivalent  result  calculated 


i 
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from  general  formulae  of  Eshelby  (1957).  If  L  Is  identified  with  one 

✓V  c 

of  the  phases,  say  L-  D  =  L  .  it  is  advantageous  to  use  (6.21)  in  the 
equivalent  form 


L  -  L  +  c 

i 


1  [(Li'ti)  fc'?l  , 


(7.25) 


whose  similarity  with  (4.24)  may  be  noted.  Clearly,  now,  the  error 


involved  in  disregarding  terms  of  order  C  is  small  so  long  as 


(L-t,)" 


-I 


is  of  the  same  order  as  ju.  ,  for  example.  Within  this  range,  fibres 
with  an  aspect  ratio  of  10  or  more  may  be  regarded  as  long.  If,  however, 

1  is  much  stiffer  than  so  that  ( L . **•  L.  ^  is  small,  terras  of 

fcx  -C  (  / 

order  £  make  a  significant  contribution  because  the  0(1)  part  of  P 
is  singular.  The  extreme  case  is  that  of  phase  2  being  rigid,  when  (7.25) 
gives 


L  *  Li  *  c*l  /c>  i 


(7.26) 


in  which  a  singular  result  is  obtained  if  the  terms  of  order  £  are 
ignored.  Explicitly, 


L  * 


(o,  0,  0 ,  f/  0 ,  o) 


(7.27) 


The  example  treated  by  Laws  and  McLaughlin  (1979)  falls  in  the  inter¬ 
mediate  range,  the  ratio  of  Young's  moduli  for  E  glass  and  their  polyester 


resin  being  around  20,  and  the  results  of  their  calculations  differed 
appreciably  from  the  infinite-fibre  predictions,  for  aspect  ratios  up  to 
about  100.  Clearly  a  similar  situation  would  occur  for  many  other 


composites . 

Real  composites  contain  imperfectly  aligned  fibres  of  variable  length. 
Laws  and  McLaughin  (1978)  made  some  approximate  allowance  for  these  fea¬ 
tures  but  no  rigorous  study  has  yet  been  made.  One  possibility  would  be 
to  obtain  both  bounds  and  self-consistent  estimates  from  (7.25),  but  with 

P  replaced  by  the  tensor  P  defined  by  (6.22).  The  difficult  task,  of 

«**  — 

course,  would  be  that  of  realistically  estimating  the  function  &  which 
defines  the  two-point  probabilities. 

7.5  A  Body  Containing  Aligned  Cracks 

A  body  containing  aligned  penny-shaped  cracks,  whose  centres  are 
distributed  uniformly  at  number  density  ,  can  reasonably  be  treated 
as  a  limiting  case  of  a  composite  with  'spheroidal'  symmetry,  the  thick¬ 
ness  to  diameter  ratio  of  the  spheroids  being  £<*<  l  •  Willis  (1980b) 
gave  the  result 

f  ~  (0,0,  0,  0,  ./*/•) 


+  JLf -  UiK*  llrf/s, 


0.28) 


so  that,  if 
stitution  of 


L  is  identified  with  L  (  L  s  O  for  the  cracks) 

-w|  *• 

(7.28)  into  (7.25)  yields  an  upper-bound  estimate  for 


L 


13* 
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Again,  the  presence  of  the  term  of  order  £  in  (7.28)  is  crucial,  since 
without  it  the  inverse  in  (7.25)  does  not  exist.  The  result  was  given 

by  Willis  (1980c).  It  is  complicated  when  expressed  in  terms  of  L.  but 

✓%/ 

takes  a  much  simpler  form  when  expressed  in  terms  of  the  inverse  M  : 


3r 


o,  i /(j*v), 

(7.29) 


where  the  relation  between  volume  concentration  and 

number  densi  tv  has  been  employed  and  the  limit  €  ~9  O  has  been  taken  with 
number  density  fixed.  Equation  (7.29)  confirms  that  the  cracks  interact 
only  with  the  stress  components  tfJ'j  ;  f  }  as  they  should.  Hoenig 

(1979)  has  performed  self-consistent  calculations  for  the  same  configura¬ 
tion.  His  estimates  agree  with  (7.29)  at  low  concentrations;  more  gen¬ 
erally,  (7.29)  bounds  the  compliances  from  below. 
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8.  WAVE  PROPAGATION 

Problems  involving  wave  propagation  necessarily  involve  spatial 
variation  and  effective  properties  cannot  be  defined  through  global 
averages,  in  the  way  outlined  for  static  problems  in  Section  2.  Instead, 
at  least  for  random  composites,  equations  which  govern  the  ensemble  aver¬ 
age  <U^  of  the  solution  to  any  particular  boundary  value  problem  are 
sought  or,  alternatively,  is  constructed  directly.  Work  published 

to  date  has  relied  on  one  of  two  methods.  The  first  involves  solving  the 
dynamical  analogue  of  (3.8)  by  perturbation  theory.  If  this  is  done  in 
the  most  elementary  manner,  it  leads  to  a  series  some  of  whose  individual 
terms  contain  integrals  that  diverge,  which  are  then  grouped  together  in 
a  way  that  eliminates  the  divergence  by  a  technique  known  as  renormaliza¬ 
tion.  This  is  a  well-known  technique  in  many  branches  of  physics;  its 
application  to  composites  has  recently  been  discussed  by  McCoy  (1979). 
Alternatively,  the  equations  can  be  solved  by  the  method  of  smoothing 
introduced  by  Karal  and  Keller  (1964) ,  in  which  an  essential  intermediate 
step  is  to  represent  in  terms  of  < *  ^  .  The  statisti- 

cally  fluctuating  term  can  then  be  eliminated  from  the  starting 

equation  to  give  an  equation  for  <»>  itself.  Karal  and  Keller  (1964) 
truncated  their  formal  iterative  solution  after  one  iteration  to  obtain 
results  valid  for  weakly  inhomogeneous  media.  McCoy  (1973),  on  the  other 
hand,  considered  the  series  to  arbitrary  order,  to  obtain  formal  expres¬ 
sions  containing  correlation  functions  of  all  orders.  Since  these  are  not 
generally  known,  however,  the  method  is  most  useful  for  weakly  inhomo¬ 
geneous  media.  An  entirely  different  method  has  been  applied  to  composites 
consisting  of  a  matrix  with  distributed  inclusions.  Each  inclusion  is 
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treated  as  a  scatterer,  which  scatters  the  field  incident  upon  it.  The 
complication  is  that  the  field  incident  upon  any  one  inclusion  has  con¬ 
tributions  from  the  scattered  fields  of  every  other  inclusion.  The  scat¬ 
tered  fields  have  always  been  represented  explicitly  as  eigenfunction 
expansions  (typically  involving  Hankel  functions) ,  so  that  the  method 
appears  to  be  restricted  to  an  isotropic  matrix,  at  least  in  the  form  pre¬ 
sented.  After  some  complex  manipulations,  equations  for  the  coefficients 
defining  the  fields  scattered  from  the  inclusions  are  derived.  These  are 
ensemble  averaged  (essentially  as  in  (6.12))  and  the  system  is  closed  by 
making  the  quasicrystalline  approximation  of  Lax  (1952)  (essentially  as 
in  (6.14)).  Particular  examples  of  this  approach  are  provided  in  the 
papers  of  Bose  and  Mai  (1973,  1974),  Datta  (1977,  1978),  Twersky  (1977), 
Varadan,  Varadan  and  Pao  (1978)  and  Varadan  and  Varadan  (1979). 

A  compressed  and  preliminary  account  will  now  be  given  of  a  varia¬ 
tional  approach.  It  is  a  natural  extension  of  that  already  given  for 
static  problems  and  offers  a  unified  treatment  for  materials  such  as 
polycrystals  (for  which  only  perturbation  theory  has  been  available)  and 
matrix- inclusion  composites  (for  which  only  the  multiple  scattering  method 
has  been  available),  without  restriction  to  an  isotropic  matrix. 

8 . 1  Variational  Formulation 

The  equation  of  motion  for  any  continuum  may  be  given,  in  the  absence 
of  body  forces,  in  the  form 

***  ~  ~  -f  )  (8.1) 

where  denotes  momentum  density  and  the  superposed  dot  represents 
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and  S.  »  M  are  the  derivatives  with  respect  to  t  of  «S  >  A7  • 

—r  t  — 

Equations  (8.9)^ and  (8. 10)^  generalize  (3.16)  and  (3.7).  The  object  now 
is  to  produce  a  variational  principle,  analogous  to  the  Hashin-Shtrikman 
principle,  that  is  equivalent  to  equations  (8.9).  One  possibility  is  to 
construct  an  analogue  of  Hamilton's  principle,  using  a  Lagrangian  that 
is  quadratic  in  £  ,  ff  .  Such  a  principle  was  displayed  by  Willis  (1980d) 
but  an  alternative,  which  appears  to  be  better,  is  to  obtain  a  principle, 
similar  to  that  of  Gurtin  (1964),  which  is  tailored  to  the  type  of  ini¬ 
tial  conditions  that  are  actually  imposed.  To  develop  this,  let  U j  , 
p  ,  ,  G.  be  displacement,  momentum,  stress  and  strain  fields  gener- 

ated  by  polarizations  with  homogeneous  initial  and  boundary 

conditions,  so  that  they  satisfy  (8.1),  (3.10),  (8.3).,  (but  not  neces¬ 
sarily  (8.2))  and  define  fields  with  a  suffix  2  similarly.  It  is  now 
easy  to  show  that  the  bilinear  form 


(8.11) 


where  it  denotes  the  ordinary  operation  of  convolution  with  respect  to 
time,  is  symmetric  and  it  follows  immediately  that  equations  (8.9)  are 
equivalent  to  the  variational  principle 


S/(t,<r)  - 
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where 

Mi,*)  -  ros-w'i-s-^J 

*■  *> [ £ -f.T V  - « -«„]/,  <8  i3> 

£  and  U  being  expressed  in  terms  of  ~C  ,  K  via  (8.4).  The  functional 
(8.13)  is  a  generalization  of  (3.17). 

The  principle  (8.12)  may  now  be  employed  to  generate  approximate 
solutions  by  the  procedure  given  in  Section  6.  If  is  replaced  by  its 

expectation  value  in  (8.12),  this  equation  generates  exactly  ’hierarchy 
equations'  such  as  (6.12),  by  suitable  choice  of  variations.  These  are 
of  no  use  as  they  stand  but  the  observation  that  they  can  be  obtained 
in  this  way  encourages  the  production  of  approximations  generated  from 
trial  fields  C  ,  n  ■  The  natural  extension  of  (6.8)  is 
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These  equations  generalize  (6.11)  and  are  the  same  as  would  be  produced 
from  the  exact  'hierarchy  equations'  by  making  the  quasicrystalline  assump¬ 
tion  (6.14)  with  regard  to  both  "T  and  If  ■  Similar  equations  were 
produced  by  Willis  (1980c),  for  a  matrix  containing  inclusions.  In 
that  work,  however,  they  were  derived  making  the  quasicrystalline  assump¬ 
tion  explicitly,  without  reference  to  a  variational  principle.  More 
general  approximate  equations  could  also  be  obtained  from  (8.12),  whose 
full  implications  are  still  under  investigation.  The  article  is  con¬ 
cluded,  however,  with  a  short  discussion  of  the  properties  of  (8.15). 

8 . 2  Wave  Speeds 

The  polarizations  (8.14)  generate  a  mean  wave 

<£>*  *  Q  £5*),  (8.16) 

obtained  by  taking  the  mean  of  (8.4).  Correspondingly,  using  (8.16), 
equations  (8.15)  can  be  recast  in  the  form 

OH.)V  -  <*>, 


(8.17) 

An  elementary  and  important  problem  for  the  ordinary  elastodvnamic  equations 


70 


is  to  determine  the  speeds  with  which  plane  waves  may  travel  in  an  in¬ 
finite,  homogeneous  medium.  The  corresponding  problem  for  (8.15)  is  to 
determine  possible  plane-wave  solutions,  for  which 


2,(z,t)  '  £r*p 


(8.18) 


with  tr ,  rr  ^  constants,  in  an  infinite,  statistically  uniform  medium. 
Thus,  plane-wave  solutions  of  (8.15)  are  required,  when  U  =  ,  the 

probabilities  are  translation-invariant  and  the  operators 

S  ,  M  take  limiting  forms  corresponding  to  infinite  V  •  The  limit 

MV 

is  easily  taken  for  the  equivalent  equations  (8.16),  (8.17).  Equation 
(8.16)  gives  (with  ia  a.  O  ) 


<u>  =  -f  0«,«)<I>  -  fiffca, «)<£>, 


(8.19) 


where 


<X>  *  2  c,  j,  s 

<£>  ’*  f  c3  25  ^f-*(ka 

S'  A/ 

and  ^  ,  PI  are  the  Fourier  transforms  of  the  operators  5  ,  Pi 

♦J  “•'*  -v' 

Equation  (8.17),  implies,  when  V  is  large 


(8.20) 
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=  -  |,da,-)f  s£  '  «/*l-“)fc.2,,(8  n) 


where  are  time-reduced  operators  obtained  from  the 

time-reduced  infinite-body  Green's  function  £  (£,«')  .  Equation 

(8.17).  can  be  expanded  similarly.  The  choice  X  —  O  could  be  made 

4  ^ 

in  (8.21)  because  of  the  translation-invariance  of  ^s|f 

Further  progress  now  depends  upon  finding  a  suitable  representation 
for  the  Green's  function  Q  .  This  satisfies  the  equatic 


ion 


which  may  be  solved  by  the  method  that  was  explained  for  the  special 
case  (4.1):  Willis  (1980b)  showed  in  detail  that 


(8.22) 
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where  Cf(i  )f  satisfy 
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(8.24) 
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so  that  they  represent  the  speeds  and  polarizations  of  plane  waves  that 


can  propagate  in  the  direction  9  in  the  comparison  material.  In 
(8.23),  the  eigenvectors  U  are  taken  as  normalized.  The  limiting 
case  CO  s  0  reproduces  the  static  formula  (4.4). 

Attention  will  now  be  restricted  to  long  waves  or  equivalently,  to 
small  to  and  a  lowest-order  approximate  solution  of  (8.21)  and  the  cor¬ 
responding  equation  for  JT  will  be  constructed  by  retaining  only  terms 
of  order  zero  in  k.  or  CO  .  Amongst  the  operators  >5  ,  M  ,  5  , 

X-  /-  K  •*-  t 

,  only  the  first  has  a  term  of  order  zero  in  tO  ,  since  the  other 

^  ^  40 

three  involve  at  least  one  time  derivative.  Furthermore,  >5  |  , 

c  * 

defined  by  (4.5).  The  Fourier  transforms  O  etc.,  on  the  other  hand, 

^  X 

are  all  homogeneous  of  degree  zero  and  are  retained  exactly.  Equation 
(8.21)  now  reduces  to 


*  f  J<sTVXP*«^ -<*)£,  “ 


(8.25) 


while  tlie  corresponding  equation  for  Tf  becomes 


(f- >)"'{[,  * 


(8.2ft) 


where  <&>  is  the  strain  associated  with  and  both  are  evaluated 
at  Equations  (8.25)  are  identical  with  (ft. 15).  Therefore,  by 
definition  of  L  ,  they  imply 


<17 


(8.27) 
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(each  side  having  the  exponential  dependence  restored).  Also,  from  equa¬ 
tions  (8.26), 


<0  -Cj -!•)<“>, 


(8.28) 


where  f  is  the  mean  density  of  the  composite  .  Equations 

(8.27),  (8.28)  show  that  the  mean  polarizations  <!>•  <E>  are 
those  that  would  be  generated  relative  to  L  ,  >  if  it16  mean 

wave  <*>  propagated  in  a  homogeneous  medium  with  moduli  1_  and  den¬ 
sity  J  '  This  conclusion  can  also  be  reached  by  algebraic  manipulation. 
Retention  of  terms  of  higher  order  in  the  equations  produces  both  dis¬ 
position  and  attenuation.  Attenuation  in  a  matrix  containing  inclusions 
has  been  studied  by  Willis  (1980c).  A  more  complete  study  of  the  pre¬ 
sent  system  is  being  conducted  and  will  be  reported  elsewhere. 
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